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ON  LIMITING  EMPIRICAL  DISTRIBUTION  FUNCTION  OF 
THE  EIGENVALUES  OF  A  MULTIVARIATE  F  MATRIX* 

Z.  D.  Bai,  Y.  Q.  Yin  and  P.  R.  Krishnaiah 


ABSTRACT 

In  this  paper,  the  authors  derived  an  explicit  expression  for  the 

limit  of  the  empirical  distribution  function  (e.d.f.)  of  a  central 

multivariate  F  matrix  when  the  number  of  variables  and  degrees  of 

freedom  both  tend  to  infinity  in  certain  fashion.  The  authors  also 

extended  the  above  result  to  the  case  when  the  underlying  distribution 

is  not  necessarily  multivariate  normal  but  the  first  four  moments 

exist.  The  limiting  distribution  is  useful  in  deriving  the  limiting 

distributions  of  certain  test  statistics  which  arise  in  multivariate 

analysis  of  variance,  canonical  correlation  analysis  and  tests  for  the 

t  f  . 

equality  of  two  covariance  matrices. 

/  .  i 

Keywords  and  Phrases:  Canonical  correlation  analysis,  empirical 
distribution  function,  large  dimensional  random  matrices,  limiting 
distribution,  multivariate  analysis  of  variance,  multivariate  F 
matrix . 


*This  paper  is  a  revised  version  of  the  CMA  Technical  Report  No.  84-42. 
In  this  version,  we  added  Section  5  which  gives  some  results  when  the 
underlying  distribution  is  not  multivariate  normal. 


1 .  INTRODUCTION 


Various  test  procedures  in  multivariate  analysis  are  based  upon 
certain  functions  of  the  eigenvalues  of  random  matrices.  A  considerable 
amount  of  work  was  done  in  the  literature  on  the  asymptotic  distribution 
theory  of  these  statistics  when  the  sample  size  is  very  large.  But, 
many  situations  arise  in  multivariate  data  analysis  when  the  number  of 
variables  and  the  sample  size  are  both  very  large.  So,  there  is  a 
great  need  to  investigate  the  distributions  of  various  functions  of  the 
eigenvalues  of  large  dimensional  random  matrices.  Distributions  of  the 
eigenvalues  of  large  dimensional  random  matrices  arise  (e.g.,  see 
Mehta  (1967))  in  nuclear  physics  also. 

Some  work  was  done  in  the  literature  on  the  limiting  empirical  dis¬ 
tribution  function  (e.d.f.)  of  large  dimensional  random  matrices.  Here, 
we  note  that  the  e.d.f.  of  a  random  matrix  Z:  p*  p  is  defined  as  N(x)/p 
where  N(x)  denotes  the  number  of  the  eigenvalues  of  Z  which  are  less  than 
or  equal  to  x.  The  e.d.f.  (also  known  as  spectral  distribution)  of  Z  is 
useful  in  deriving  the  distributions  of  certain  functions  of  the  eigenvalues 
of  Z. 

Now,  let  Sj  :  p  *  p  be  distributed  as  central  Wishart  matrix  with  m  degrees 

of  t reedom  and  E(S^/m)  =1^.  Also,  let  p  and  hi  both  tend  to  infinity  such  that 

lim(p/m)  ■  y  >  0.  Then,  it  Is  known  (see  Crenander  and  Silversteln  (1977), 

Jonsson  (1982)  and  Wachter  (1978))  that  the  e.d.f.  F  (x)  of  the  eigen- 

m 

values  of  S./m  tends  to  F  (x)  where  F  (x)  is  the  distribution  function  with 
1  y  y 

density  function  given  by 


f  (x) 

y 


a)  (b 
2irxy 


x) 


a  <  x  <  b 

otherwise 


2  2 

where  a  -  (l  -  /y)  ,  and  b  -  (1  +  /y)  and  0  <  y  _<  1;  for  I  <  y  <  »,  F^(x)  has 
mass  1  -  (1/y)  at  zero  and  f^(x)  on  (a,b).  Yin  and  Krishnaiah  (1983b)  showed 
that  the  spectral  distribution  of  the  sample  covariance  matrix  has  a  limit  when 
the  underlying  distribution  is  isotropic  and  y  <  1. 


Yin  and  Krishnaiah  (1983a)  showed  chat  the  spectral  distribution  of  S ^ T /  m 

tends  to  a  limit  in  probability  for  each  x  under  the  following  conditions: 

(a)  T  is  a  symmetric,  positive  definite  matrix  and  G  (x)  is  the  e.d.f. 
of  the  eigenvalues  of  T, 


(b) 

(c) 

(d) 


and  T  are  independent  of  each  other, 

lim(p/m)  =  y  exists  and  finite 
p  ,nr*-°° 

‘x^dG  (x) — >-8.  exists  in  for  k  =  1,  2,...  and  EH-^ 
J  p  k  4k 


Yin  and  Krishnaiah  (1984b)  extended  the  above  result  to  the  case  when  is  the 
sample  sums  of  squares  and  cross  products  matrix  based  upon  observations  from  an 
isotropic  population. 


Now,  let  S^:  px  p  and  S„:  p*  p  be  distributed  independently  as  central  Wishart  matr 

with  m  and  n  degrees  of  freedom  and  K(Sj/m)  ■  lifS^/n)  “  Ip.  Then,  the  distribu¬ 
tion  of  nSjS'Vm  is  known  to  be  the  central  multivariate  F  matrix.  Applying  the 
result  of  Yin  and  Krishnaiah  (1983),  Yin,  Bai ,  and  Krishnaiah  (1983)  showed  that  the 

limit  of  the  spectral  distribution  of  the  central  multivariate  F  matrix  exists 
when  p/m-^+y  exists  and  (p/n)  — *  z<  h  as  p  -*■  “ .  Silverstein  (1984b)  showed  the 
validity  of  the  above  result  even  for  the  case  ^  <  z  <  1  by  making  a  minor 
modification  in  the  proof  of  Yin,  Bai,  and  Krishnaiah  (1983).  Yin  and  Krishnaiah 
(1983a) gave  an  expression  for  the  moments  of  the  limit  of  the  e.d.f.  of  the 
eigenvalues  of  S^T/m.  Starting  from  this  expression,  Silverstein  (1984^)  derived 
an  explicit  expression  for  the  limit  of  the  e.d.f.  of  the  eigenvalues  of  the  multivariate 
F  matrix.  Wachter  (1980)  had  earlier  obtained  the  explicit  expression  by  using  a 
different  method.  In  this  paper,  the  authors  give  an  alternative  derivation  of 
the  above  method.  The  authors  also  gave  explicit  expressions  for  the  moments  of 
the  above  limiting  distribution  and  these  expressions  are  not  known  in  the 
literature.  Finally,  the  authors  proved  the  existence  of  the  limit  of  the  spectral 
distribution  of  nS^S;  Vmwhen  the  underlying  distribution  is  not  multivariate 


normal  but  the  first  four  moments  exist. 
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2.  PRELIMINARIES 


In  this  section,  we  give  some  results  which  are  needed  in  the  sequel  as 
well  as  a  brief  review  of  known  results  on  limiting  spectral  distribution 


of  a  multivariate  F  matrix. 


Lemma  2.1  Let  z 


then 


e  (.0,1)  ,  a'  =  (1  -  and  b'  -  (1  +  /z)2.  If  0  j,  I  cl  <  a'  » 


J_  b  -7 _ 1 - v  {(x  -  a’)(b'  -  x)P  dx  -  y-[l  +  z  -  t  -  {(1  -  z  -  t)2  -  4tz}^ 

2nz  J  ,  (x  —  t)  zz 

Ja  (2.1) 


Proo l  Making  the  transformation  u  =  [ 2x  ~(b'  +  a')]/(b'  -  a')  in  the  left 
side  ot  12.  n  ,  we  obtain 


K(t)  = 


2nz 


rb'  l  A  (b'~  a)  fl 

tTri-rT((x-a'Xb'-x))V--^ 


r!  / 1  2A 

(1  -  -U.J-  du 

J_!  u  +  A 


(2.2) 


wiie re  Z  =  (b'  t-  a1  -  2t)/(b'  -  a').  It  is  known  (see  Jonsson  (1982)) 
that  R ( 0)  =  1.  So,  for  any  r  r  (0,1),  we  have 


i  /  r 


1  , ,  2 .  Hi  , 

_ (1  -  u  )  du 

-1  u  +  ((1  +  r)/2/r  ) 


1. 


(2.3) 


Now ,  let  &  =  (I  +  r ) / 2*  r .  Since  A  >  1,  the  condition  r  c  (0,1)  is  satisfied, 
So,  using  (2.2)  and  (2.3),  we  obtain 


3  2  U 

A“  -  l)/4z  -  [1  +  z  -  t  -  (d  -  z  -  t)  -  4tz }  ]/2z. 

(2.4) 

Lemma  2 . 2  For  any  nonnegative  integers  m  and  w,  we  have 


[w/2] 


(-1) 


J  =0 


2m  -  2  ?^  0w  f  m 
2m  -  w 


(2.3) 


Here 


defined  to  be  zero  when  m  <  w. 


P.root~  If  m  ■  0,  the  proof  is  trivial.  We  now  prove  the  result  for  m  >  0  by 
induction.  Suppose  the  result  is  true  for  a  fixed  value  of  m.  Then,  we  have 


[w/2] 

l 


2=0 


(-1) 


/  2m  +  2  -  2?\ 

\^2m  -  w  +  2  J 


m+1 


l  (-D 

lmQ 


l 


m  ^  /2m  +  2  -  2  r  j 

-  \j  \^2m  +  2  -  w  J 


-  j  (-nl 

t=0 


C:')- 


So  the  result  follows. 
Lemma  3 


(ttY  +  fl)dy 

J  (v  +  cy  -  oy2) 3^2 


=  4QY  ~  2B6  +  (AqB  +  2o_6)y 

(^^+4yc)(y  +  'y-f.  y2) 


+  const. 


(2. -6) 


The  above  lemma  can  be  verified  directly  by  differentiation. 
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3.  MOMENTS  OF  THE  LIMITINC  SPECTRAL  DISTRIBUTION 
OF  A  MULTIVARIATE  F  MATRIX 


Let  F  (x)  denote  the  limit  of  the  e.d.f.  of  the  multivariate  F  matrix 

y.z 

CO 

as  defined  in  the  preceding  section.  Also,  let  {E^}  ^^denote  the  moments  of 
F  (x).  Then,  from  Yin  and  Krishnaiah  (1983),  we  know  that 

y.z 


E*  ‘  J/CO'dT)' 


(3.1) 


where 


B(k«w)  “It-? 


(w  +  1) !  „nl  A-w 
n^  ...  n^!  1  *•*  k-w 


(3.2) 


and  the  summation  in  (3.2)  is  over  all  possible  values  of  nj,...,n^_w  subject 

to  the  restrictions  n.+  •••  +n,  =  w  +  1  and  n,  +  n~  +  •••  +(k  -  w)rL  =  k. 

1  k-w  1  2  k-w 

Also,  =  E ( x” 1 )  for  i  =  1,  2,...  where  the  density  of  x  is  given  by 


gz(x)  = 


(x  -  a')(b’  -  x)  a,  <  x  <  b. 
2irxz 

)  otherwise 


(3.3) 


where  a'  ,  b'  ,  and  z  are  as  defined  in  the  preceding  section.  For  any  |t|  <  a'  , 


2  ') 

4>(t)  =  tHj  +  t  ll2  +  t  ll.3  +  .  .  .  . 

=  E(t/ (x  -  t)  ) 

1  fb' 

=  ~  ( t  { (x  -  a '  )  ( b '  -  x)}‘/U  -  t)]dx 

2ttz  , 
a 

=  (1  -  z  -  t )  f  1  -  (1  -  { (4tz/  ( 1  -  z  -  t)2)}ii)]/2z  (3.4) 

by  using  Lemma  2.1.  Since  B(k,w)  is  the  coefficient  of  t  in  Taylor  expansion 
of  {y(t))W+^,  we  obtain 


(3.5) 
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+  /  0  implies  0  ^  j  +  1  -  k_<0,  l.e.  j  «  0. 

—  '  k  -  j  j 


w-t-l  /  ,\ 

Noticing  V  I  1  .-l)1  *  o,  we  see  that  the  expression  of  B(k,w)  only  contains 
<.=0  \  ’  / 

f  7  /t\ 

the  terms  with  j  >  <(w+l)/2).  If  l  is  even,  then  j  /  =  ®  implies 
j  •-  !2  (w+l)/2,  which  is  contrary  to  j  >  ((w+l)/2).  Thus,  in  the  expression  < 

B(k,w),  tliere  are  only  the  terms  with  j  >  ((v+l)/2)  and  being  odd.  Applying 
Lemma  2.2,  we  obtain 


ls(k,w)  =  £ 


iv^21  /  w  +  1 V  i  Nk+1  - 

U-l  ) 


\ 


2r+l. 


?  +  1/2  j  22j-l-wzj-w-l 
j=( (w+3)/2)  V.  J 


X  ,rw+1:2^  (1  -  z)^l“j"k 

k  ~  J 


l-«]/u+1N  f  (2t  +  1)  ! (2.1  -2-201 


1  (“-■ 

im 


0  V:+lJ  J  =  ((w 


(  (w+3) /2) 


j  \l  1  ( j  -  l  -  1)  ! 


(-nVV-^u  -  (*  +  j  :  «  - 


V  (k  +  j  -  w  -  2) ! (w  +  1)  j -w- 1  . .  .vrt-l-j-k 
1  -i!(j  -  «  *  l)!(k  -  j)*  Z  U  *  z) 


j=w+l 


k  y  1  _ -i-"  J2 _  ZJ(!  _  z\ 

J (|+w+!)!j!(k  -  j-w-l)r  U  z> 

j«=U  J  J 


-k-j 


(3.6) 


•rum  0.j)jnd  0.6)  it  follows  chat 
k-  1  /,  \  k-w-  1 

i 
i 


l., 


k- 1  /.  \  k-w- 1  , 

JoV  (w)  1  ~jY(j  +  w  + 


(k  +  j  -  1)! 


1) f(k  -  j  -  w  -  I)f 


J  (1  -  z)  k  ^ 


k-I-j  , . 

»  l  *Jo  -  <rk‘J  i  v(k 


k-  l 

l  > 

j-U 


(k  +  j  -  1) 


w=Q 


wy  j !  ( j  ♦  w  +  l)!(k  -  J  -  w  -  1)!  * 

(3.7) 


The  k-th  moment  o I  I  (•)  can  be  easily  computed  from  (3.7) 

y»z 


7 


4.  DERIVATION  OK  THE  LIMITING  SPECTRAL  DISTRIBUTION 
OK  MULTIVARIATE  F  MATRIX 


Substituting  the  well-known  formulas 


_  V  ^k-j-w-l'\/j+w+l\ 

t=0  V  c  /  V  w  "  c  . 


(4.1) 


into  (3.7),  we  obtain  the  following  by  changing  the  order  of  summation: 

k;1  i  _k  ,  [(k-1-3 )/2]  .  v  ^  \ 

Ek  =  l  zJ(l  -  z)  k_J  l  (k  +  j  -  1) !  y*  (1  +  y)k~3  1-2t 

j=0  t=0  (k  -  .1  -  1  -  2t)  !  t !  (j  +  t  +  l)!j! 


k-1  . 


‘V*  i  _k_i  Kk-j-l)/2]  t-i  k+-s_l_2t 

=  [zJ(l-z)kJ  Y  (k  +  j  -  1)  !y  J(I  +  y)**3 

j=0  t=0  (k  +  j  -  1  -  2 1 ) ! ( t  -  j) ! j ! (t  +  1)! 

.  Y  (2t\  _l  J  zi(I  .  z)-k-j  /cYk  +  j  -l'.  yt-J(l  +  ,k+3-l-2t 
t‘o  V  /  c  +  1  j-o  VLV  2t  . 


Using  the  formulas 


in  (4,2),  we  obtain 


K  t=0  \  t  /  t  +  i  s=c  ^  s  Jj=o 


(4.2) 


(4.3) 


t-j..  .  >.  k+j-  1  -  2 1 

y  (1  +  y) 


t  “  Y  Vs -t 


W  A  V 


s:Y  z^d  -  z)'k-t+j 


j  ,  ,  k—  t—  1  —  j 

r  (l  +  y)  J 


2t-s,.  N-k-t 
z  (1  -  z) 


/ 1  ,  .  k  s  1  z  v S—  C 

(1  +  y)  (y  +  z) 


X 


k~Y2t\_j_  v  i k_  1 Y c  ^  zc~ 

j-A  +  \  _LU + 1  { s 


s=0  'v 


-s,.  ,-k-t 

( 1  -  z) 


,,  ^  , k-t-s-1  ,  ,s 

X  (1  +  y)  (y  +  z) 


(A. 4) 


■ T2 '  T(ks;li Ok'' 


s  =  0  t=s 


,  ,  k-  s  -1-t  .  ,  ,  s 

X  (1  +  y)  (y  +  z) 


=  l(k-y/2J  k  \  29  fk-l\fk-l-2s\  (2s)  !  (2s  +  2t)  ! 

s=0  c  =  0  \  2s  \  C  ,  s!(2s  +  t)!(s  +  c  +  1) 


t,.  .-k-t-s.,  ,  ,k-2s-l-t.  ,  .s 

*  z  (1  -  z)  (1  +  y)  (y  +  z) 


(4.5) 


2  2  —  U 

Define  a  random  vector  (U  V)  where  the  marginal  density  of  U  is  —  (1  -  x"")  ‘I 

’  TT  I 

and  the  conditional  density  of  V  given  U  =  x  is 


— i-L2  ,,  ->J/2  ...  0,1/2  (v). 

1  -  l ( 1 -X- )  , ((1-X-)  ] 


(4.6) 


It  is  oasv  to  see  that 


F  .  2s  2t+2s 
E  t  b  V 


( 2s) ! (2t  +  2s) ! 
s!(2s  +  t) ! (s  +  t  +  1)! 


(4.7) 


and  that 


t:{u2s+tv2t+2s+1 )  =  0. 


Hence,  t rom  (4.5),  (4.7)  it  follows  that 


k- 1  k-l-s 


Eu  =  E  l  l 


j=0  t=0 


k  -  1  k-l-s  t 


z  (1  -  z) 


,.  .  xk-s-t-1,  .  ,2  ,,s„2t+s 

x  (1  +  y)  (y  +  z)  U  V 


(1  -  z) 


s  11  + > +  r^E + 


4V2z  4UV/y  +  z  \k_1 


(4.8) 


Now,  we  compute  the  distribution  of 
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,  ,  .  4V  z  ,  4 UV  / y  +  z  1  ,  „ 

1  +  y  +  i +  - - -  /  (1  -  z)  . 

1  '  Z  /T^T  . 


(4.9) 


Let  Wj  =  UV  and 
of  (w^  ,  w.-,)  is 


=  V  .  Then  we  can  easily  show  that  the  joint  distribution 


1  t  2, -3/2 

-  w0(w_  -  w.)  I  9  9 


n  2  2  1- 


[w“+w0  w2  . 


where  I  [ y  ]  ta^-es  value  1  if  o  <  B  is  true  and  zero  otherwise.  Applying 

_  .1/2  ]  >2 

Lemma  2.3,  we  can  compute  the  density  of  4[-  (y+z)  *w/(l-z)  '  +  z  w,/  (l-z  i' ] 

as  follows: 


Let  q  <  q.,  be  the  roots  of  the  equation 

(1  -  z)2  2  ,  1  - 


q  +Tb-(7TTT  (x"  q) 


2  1  -  z 


(4.10) 


in  the  variable  q.  Let  a  =  1  ,  6  -  0,  .  =  {(1  -  z)/16(y  +  z)}x  ,  <S  *  -'((l-z)/4z)  * 
(( 1  -  z)x/8(y  +  z) ) :  and  .  =  ( 1  -  z) / 1 6 ( y  +  z) .  Then 

Al  "  4“'  -  2BS  •  -  4  (y  1 


A 2  =  4  Bp  +  2nd 


I  -  z  1  -  z 

2z  4(y  +  z)  x 


A3  =  S2  +  4yp  -  li-  .-2.)/[y  +  z  +  x-z1  . 

lbz  (y  +  z) 

1/2  1/2,  , 

By  Lemma  2.3,  we  get  the  density  of'Aw^ty+z)  “"/(l-z)  “  ■  + ' 4zw0/(l-z) •  is  as 

given  below: 


f,(x)  .  I -Z3Z3  u_i*2l  f 

4/y  +  z  16z^  J c 


(  1  -  z  1  -  z  ,  .2 

q  (  4z  q  1 6 (y  +  z)  X  q 


0  -  z )• 


Aj  +  A2q 


64nz 2/y  +  z  A,/(l  -  z)  1  -  z  7  72  q, 

7~<r~  q  ~  i6(7T'ou~  ^  1 

r  n  proving  (4.11),  we  use  the  formula  fy+^(x)  =  j  f (x-q ,q)  dq  where  fx+^(.)  and  1  (. 
are  the  densities  of  x+y  and  (x,v)  respectively.  Note  that  the  integrand  in  the 
integral  in  (4.11)  is  zero  outside  the  interval  [  q  ^ ,  q  ^  ]  by  the  indicator  factor. 
Since  q  and  q.,  are  roots  of  equation  (4.10),  we  have 


,  -  ^  f - .  _  \  1  i-  ) 

64ttz  /y  +  z  j 


:i  -  z)5/2  ~,Z  (  V  (q2  ~  qi)  (4.12) 


2 y— TT  a-,  (1-2)  q,q-, 


64itz  / y  +  z  "3  '  MP2 

Here  we  note  that  42  q  -16TF?)  <!“q)2  ■  <  qr  £or  q  ■  qi or  q: 


From  (4.10),  we  can  compute 


1  -  y  +  z  -  yz 


(4.13) 


2zJ(y  4-  z)(-x  (1  -  z)  4-  4X2  +  4(y  4-  2)) 
'l  y  +  z  -  yz 


(4.14) 


From  (4.12)  -  (4.14),  it  follows  that 

5/9  4  [x2 (1-z) /  4 (y+z)  ]  /(y  +  z)  (-x2  ( 1  -  z)  +  4xz  +  4(y  +  z)) 

,  ,  (1-z)  “  12Hz  (y-t-zl _ » _ 

t,(x)=- - y  ■  "  -  .2  0  9 

1  64*zVv  +  z  <1  -  *>  [y  +  2  +  xz](l  -  z)  z‘r 


r  —  -  - - - -  ■  -  ~  - 

————7 - r  /-x“(l  -  z)  +  4x2  +  4  (y  +  2)  . 

1  +  z  +  xz) 


(4.15) 


From  Lhi.s  we  can  easily  obtain  the  density  of  1  +  y  +  -j- - - 


4UV/7+  Z 


l  ,(x) 


(I  -  z)f-(l  -  z)~x"  +  2(1  -t  v  -F  z  -  vz)x 

2n(xz  +  y) 


-  (1  -  vV 


(1  -  z ) ^  (x  -  a) (h  -  x) 
2^(xz  +  y) 


(4.16) 


( 1  -  /v  +  z  -  yz  f  ,  , 
where  a  - - 1 L -  and  b 

(1  -  z) 


(1  +  /y  t  z  -  yz  f 
(1  -  z)2 


Since  f^(x)  ^  0  if  and  only  if  equation  (4.10)  has  two  different  roots, 
we  find  by  checking  the  steps  of  computation  that  f2(x)  4  0  if  and  only  if 


a  <  x  <  b.  Recalling  (4.9),  we  obtain 


Ek  "  TTT  j  xk-lf2(x)dx 
J  a 

fb  k[  f2(X>  1 

“  L  x  -  *>J 


(4.17) 
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Now,  let 


f(x)  = 


f (x) 
x(l  -  z) 


(1  -  z)/ (x  -  a)  (b  -  x)  . t_ 

- - 2nx(y  +  xz) - If  a  <  X  <  b, 

0  otherwise. 


(4.18) 


In  (8],  it  is  shown  that  the  distribution  F  _  is  determined  by  all  its  moments. 

y  *  z 

From  (4.11)^4.12)  it  follows  that 


Fy.Z(x)  "  Yl{(0,»)}(x)  +  J_  f(x)dx’ 


(4.19) 


where  y  »  1  -  f (x)dx 


Finally,  we  only  need  to  compute  the  integral  f(x)dx.  If  we  set 


u  =  [2X  -  (b  +  a) 1 / (b  -  a),  we  get 


I  =  f(x)dx  = 


(1  -  z) 


/(x  -  a)  (b  -  x)  ,  i  z  /(x  -  a)  (b  -  x)  , 
- dx  -  - - 7 — ; - 7 - rr - dx 

x  j  (z(x  +  v) ) 

a  Ja 


»  -  zHb_-_a)  [f1  A  -  J*?  du  .  [b  /  du 

47ry  J_j  u  +  X1  Ja  u  +  X2  J 


,  ,  b+a  ,  .  b  +  a  +  (2y)/z  „  .  t. 

where  X  =»  -r -  and  X-  *  - .  —  —  ^ - .  Using  (2.3)  we  get 

1  b  -  a  2  b-a 


T  “  (I  ~  ZUb'  -"  <X!  *  AJ  -  1)  -  (X2  -  /x| 


(1  -  z) 


2  1/2  -r - 


[(z^ab  +  zy(a+b)  +  b  )  /.z  -  vab  — (y / z) ] 


-  z)  y  +  z _ 1 1  -  y  1  _  v 

2y  z(l  -  z)  1  -  z  z 


-  -J7  [1  +  y  -  U  -  yll 


if  0  <  y  <  1. 


if  y  >  1 


Hence 


if  0  <  y  £  1 . 


1  -  -  if  y  >  1. 

y 


Substituting  this  into  (4.19),  we  get  the  expression  of  F  z(x) 
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5.  LIMITING  DISTRIBUTION  WHEN  THE  UNDERLYING  DISTRIBUTION  IS  NOT 

MULTIVARIATE  NORMAL 

In  this  section,  we  prove  the  existence  of  the  limiting  spectral 
distribution  of  the  random  matrix  F  when  the  underlying  distribution  is 
not  multivariate  normal  but  the  first  four  moments  exist.  We  need  the 
following  difinition  and  lemmas  in  the  sequel. 

Definition  5.1.  A  random  vector  x'  =  (x. , . . . ,x  )  is  said  to  be  M  -  PD 
-  1  P 

(projection  distribution)  bounded  if  there  is  a  positive  constant  M 

such  that  for  anv  constant  unit  vector  a'  =  (a  ,...,a  )  and  any  x  >  0 

1  p 

Pt'a'xlj^xl^Mx  (5.1) 

Definition  5.2.  A  sequence  of  random  vectors  x^  of  order  p(n)  *  1  is 

said  to  be  M  -  PD  bounded  if  there  is  a  common  constant  M  such  that  each 

vector  x  is  M  -  PI)  bounded. 

~n 

Example  5.1.  If  x  is  distributed  as  a  multivariate  normal  with  mean 

(Y 

vector  0  and  covariance  matrix  I,  then  x  is  y-  -  PD  bounded. 

Example  5.2.  Let  x  be  uniformly  distributed  on  the  p-dimensional  sphere 
with  center  at  zero  and  radius  >p,  p  jl  3 .  Then  x  is  2  -  PD  bounded. 

Lemma  5,1.  Let  x  be  an  M  -  PD  bounded  random  vector  and  let  y  be  a 
random  variable  with  E(l/;yj)  <  m.  In  addition,  we  assume  that  y  and  x 
are  distributed  independent  of  each  other.  Then  v  x  is  ME ( 1 / ! y j )  -  PD 
bounded . 

The  proof  of  the  above  lemma  follows  immediately. 

Lemma  5.2.  If  A:  p *  p  is  any  nonnegative  definite  matrix  and  ■  is  the  largest 
eigenvalue  of  A,  then 

|x'A  x  -  y  'Ay  !  £  >•  II  x-y  ||  2  +  21 1/2  ||  x-y  ||  Ily'Ay]!1^2  (5.2) 

where  ||x||  denotes  the  Eucledian  norm  of  x. 

The  above  lemma  can  be  proved  by  applying  Schwarz's  inequality. 
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Lemma  5.3.  If  0  <  r  <  1/2,  then  the  unit  ball  in  Rp  can  be  covered  by 

balls  of  radius  r  in  such  a  way  that  the  number  of  these  smaller  balls 

? 

does  not  exceed  C  exp{(p/2)  log(87:e/r“) }  where  C  is  a  constant. 

The  proof  of  the  above  lemma  is  given  in  Yin,  Bai  and  Krishnaiah 
(1983)  . 

Lemma  5.4.  Let  (x^, . . . ,x  , . , be  a  sequence  of  random  vectors  which  are 

M  -  PD  bounded  with  a  common  M  and  x^  is  of  order  p *  1 ,  Also,  let 

x  , .  . .  ,x  be  a  sample  of  n  =  n(p)  independent  observations  on  x  . 

-Pi  ~Pn  ~p 

Then,  for  any  a;  p*  1  and  e  >  0 


P  [  y  (a  'x  . )  ^  <_  ne  ]  <_  C  exp(  y  log  (weM^  e/2); 
i-1  -  ~pl 


(5.3) 


Proof.  Let  Y.  =  a'x  .  for  i  =  l,2,...,n.  Then  Y, , . . .  ,Y  are  i.i.d. 
-  i  «  „pi  In 

random  variables  and 


But,  we  know  that 


V(x)  =  P  [  ]  Y x  j  <_  x]  <_  Mx 


[n/2  ] 


(5.4) 


P[Yj  +  ...+Y^  <  nE]  <  P[  l  ^Y2k-1  +  Y2k}  -  ne]  (5.5) 


where  [a]  is  an  integral  part  of  a.  We  have 
P[Y^+Y^  <_  x]  =  f  V(/x^y)dV(/y) 

1  z  0  rx 

M  *  x-y  dV(v  y) 

2  °,X 

4~  f  ^y/  (x-v)  dy. 


(5.6) 


Cy/Gcy))1^2  dy  =  -~- 


(5.7) 


P [Y^  +  <  x]  <  ~ 


x . 


(5.8) 
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Now,  let  Zk  =  Y^k_1  +  Y”k .  Then  Z^. 
random  variables  a  =  [n/2]  and 


P[Z  i  x]  <  M  x 


where  M*  =  t M-/  4  .  Now  let 


F  (x)  =  P  [Z,  +  .  . .  +  Z  <  x] 
a  1  a  — 


We  shall  now  prove,  by  induction,  that 


F  (x)  <  Maxa/a!  . 
a  — 


,Z  are  i.i.d.  nonnegative 

a 


(5.9) 


(5.10) 


(5.11) 


The  above  result  is  true  for  a  =  1.  Here  and  in  (5.12)  below,  we  use 

M  instead  of  M*  for  simplicity.  We  will  assume  the  result  to  be  true  for  a =  b. 

rX 

Now  F.  . ,  (x)  =  |  F  (x-y)dF,  (y) 


£  M  ,  (x-y)dF  (y) 
J0 


=  M  ]  F  (v)dy 

J0 


rfb+i  rx 


y  dy  = 


Mb+1xb+l 
(b+1) ! 


(5.12) 


So,  by  induction,  the  inequality  (5.11)  holds  good.  Now,  using  in¬ 


equalities  (5.5),  (5.8)  and  (5.11),  we  obtain 

P[Y"r+...+Y2  <  ne]  <  (it  M2/4)a  (ne)a/a!  . 
i  n  — 


(5.13) 


Now,  by  applying  Stirling's  formula,  we  obtain  the  desired  result. 

We  now  prove  the  following  result 

Theorem  5.1.  Let  ,  p  =  1 ,2 , . . .  be  M  -  PD  bounded  and  (p/n)  -*•  y  e  (0,1). 
Then,  for  any  R  >  0,  there  exist  positive  constants  C  <  00 ,  D  <  00 , 
a  >  0,  tp  >  0  depending  only  on  y,  M  and  R,  such  that 

Pie  i  ip  1  1  R]  I  CDPeap,  0  <  e  <  eQ 

where  X  and  X  are  the  smallest  and  largest  eigenvalues  of  W 
~p  P  P 
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respectively  where 


1  “ 

W  =  -  X  .  X '  . 

P  n  -PJ-PJ 


(5.14) 


and  X  , , . . . ,X  are  independent  observations  on  X  . 

-Pi  -pn  K  -P 

Proof.  Let  S  denote  the  unit  sphere  in  R  .  Then, 

-  p  p 

a.  =  min  y'  W  y.  (5.15) 

P  v  e  S  ~  p  - 
P 

Let  r  <  1/2  be  a  positive  constant  to  be  chosen  later.  According  to 

Lemma  5.3,  let  B  (x,,r),...,B  (x  ,r)  be  those  balls  with  radius  r  and 
p  1  p  q 

centers  x^,...,x^  which  cover  the  unit  sphere  S  .  In  addition,  q  satisfies 


q  <_  C  exp{-|p  log 


(5.16) 


Then  we  have 

-9 

P(2_  1.  e  >  R)  <  /  P(  min  y '  W  y  £  e  ,  x  <  R)  . 

P  P  k=l  ytB (x^, r)0  3^  P  P 

If  y  e  B(x.  ,r)f)S  and  y'W  y  <  e,  X  <  R,  then  by  Lemma  5.2,  we  have 
-  k  p  i  p  1  —  p  — 

x,'  W  x,  <  y '  W  y  +  J  x,'  W  x,  -y'W  y ! 

.k  p  „k  —  L  pi  _k  p  „k  i  P  -  ' 


<  e  +  R  r2  +  2R1/2re1/2  -  (E1/2  +  R1/2  r>2. 


Therefore , 


P(X  <  e,  X  <  R)  <  }  P(x W  x,  <  (t1/2  +  R1/2  r)2) 

— p  —  p  —  —  _  k  p  ,k  — 


But  by  Lemma  5.4,  we  have 


P(x'  W  x  <  (e1/2  +  R1/2  r)2) 
^  k  p  „  k  — 


,  r  f n  ,  ,  re  2.  1/2  1/2  .2,,.  .2., 

£  C  exp{  —  log (  —  M  (c  +  R  r)  /(1-r)  ) 


If  we  take  r  =  R  ^2e^2,  we  have  from  (5.16)  and  (5.17), 


(5.17) 


PQ.  <  e,  X  <  R) 

p  -  p  - 

£  C  exp{^-  log  6^-  +  log(8ire  M^e)} 

<  C  £ap  DP 

for  p  large  enough  where  D  =  (8ireR)^^(8iTM^)^^,  0  <  a  <  1/ 2 ( 1/y  -  1). 

-1/21/2  A 

But  (5.17)  holds  provided  that  r  =  R  z  <  1/2,  i.e.  z<z^='z  R. 

So,  the  theorem  is  proved. 

Theorem  5.2.  Suppose  that  the  assumptions  of  Theorem  5.1  hold.  Also, 

we  assume  that  X'.  =  (X. X.  ),  and  X..,  i,j  =  1,2,...,  are  i.i.d. 
-pi  il  ip  ij 

2 

random  variables  with  common  mean  zero  and  variance  a  and  finite 
fourth  moment.  Then  there  exist  a  positive  number  e  such  that 

lim  A.  >  e ,  a.s. 
p-x»  c 

00 

2  _ 

Proof .  Take  R  >  (l+/y)  a  and  set  A  =  (J  (X  R)  .  According  to 

P  m=p 

Yin,  Bai  and  Krishnaiah  (1984),  P(A^)  -*■  0  as  p  ->  00 .  Thus  for  any 
e  <  Eq,  we  have 

00 

P(U  (A  <  e))  <  y  P(i  <  e,  r  <  R)  +  p(A  ) 

~m  —  L  m  —  m  —  p 

m>_p  m=p 

00 

<_  C  l  (eaD)P  +  P (A  )  0,  p  ->■  “>. 

m=p  P 

We  will  now  prove  that  the  limiting  spectral  distribution  of 

F  matrix  exist  under  conditions  weaker  than  assumed  in  Sections  3  and  4 

Theorem  5.3.  Let  W  =  -  X  X\  W  =  -  Y  Y'  where  X  =  (X.  .):  p*  n, 
-  pnpppmpp  p  ij 

Y  =  (Y .  , )  :  p  x  m  where  n  =  n(p),  m  =  m(p)  such  that  (p/n)  -*•  v  z  (0,™)  , 

P  iJ 

(p/m)  y'  z  (0,1).  Also,  we  assume  that  X^  ,  Y^t  (i  =  l,2,...,p; 
j  =  l,2,...,n;  t  =  l,2,...,m)  are  i.i.d.  random  variables  with  common 
mean  zero,  variance  one  and  finite  fourth  moment.  Also,  let 


F  =  W  W  where  W  is  nonsingular.  Also,  for  each  p,  we  assume  that 
P 
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(Y^ ,  .  . .  , Y  is  M  -  PD  bounded.  Then,  as  p  ■»  * ,  has  a  limiting 
spectral  distribution  with  probability  one. 

Proof .  Applying  Theorem  5.2,  we  have 

0  <  e  <  X  <  X  <R=(1  +  v7)2  +  1 
-  -p  -  P  “ 

for  sufficiently  large  P,  with  probability  one.  Thus,  for  any  k  _>  1 , 

1  k  f°°  -k  f  ^  _k 

-  tr(W  )  «  x  dQ  (x)  =  I  x  dQ  (x) 

p  p  Jo  ?  4 

where  Q  is  the  e.d.f.  of  W  . 

P  P 

According  to  Yin  (1984)  ,  the  random  matrix  W^T  has  a  limiting 

spectral  distribution  if  W  satisfies  the  condition  of  the  theorem  and 

P 

T  satisfies  the  conditions  (a) ,  (b) ,  (c)  stated  on  page  2  and  the 
following  condition: 

(d')  [  x^dq  (x)  H.  a.s.  for  k  =  1,2,...  and  Z  H  =  00 . 

j  p  k  ZK. 


Taking  T  -  I,  the  identity,  we  know  that  Qp  Q,  a.s.  where  Q  is 
the  distribution  whose  density  is  given  by  (3.3).  Thus 

rb' 


-  tr(W  )k 
P  P 


Hk  = 


x  g  (x)dx,  a.s, 
z 


Again,  using  Yin  (1984),  with  T  =  (W^)  \  the  theorem  is  proved.  In 
view  of  Yin  0-984)  and  Yin  and  Krishnaiah  (1983a)  and  the  derivation  in 
54,  we  know  that  the  limiting  spectral  distribution  is  the  same  as 
that  given  in  54. 
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